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Abstract 

Acoustic black holes are fluid dynamic analogs of general relativis- 
tic black holes, wherein the behaviour of sound waves in a moving 
fluid acts as an analog for scalar fields propagating in a gravitational 
background. Acoustic horizons, which are intimately related to re- 
gions where the speed of the fluid flow exceeds the local speed of 
sound, possess many of the properties more normally associated with 
the event horizons of general relativity, up to and including Hawking 
radiation. Acoustic black holes have received much attention because 
it would seem to be much easier to experimentally create an acoustic 
horizon than to create an event horizon. In this note we wish to point 
out some potential difficulties (and opportunities) in actually setting 
up an experiment that possesses an acoustic horizon. We show that in 
zero-viscosity, stationary fluid flow with generic boundary conditions, 
the creation of an acoustic horizon is accompanied by a formally infi- 
nite "surface gravity", and a formally infinite Hawking flux. Only by 
applying a suitable non-constant external body force, and for very spe- 
cific boundary conditions on the flow, can these quantities be kept fi- 
nite. This problem is ameliorated in more realistic models of the fliiid. 
For instance, adding viscosity always makes the Hawking flux finite 
(and typically large), but doing so greatly complicates the behaviour 
of the acoustic radiation — viscosity is tantamount to explicitly break- 
ing "acoustic Lorentz invariance" . Thus, this issue represents both a 
difficulty and an opportunity — acoustic horizons may be somewhat 
more difficult to form than naively envisaged, but if formed, they may 
be much easier to detect than one would at first suppose. 
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1 Introduction 

Acoustic black holes are very useful toy models that share many of the funda- 
mental properties of the black holes of general relativity, while having a very 
clear and clean physical interpretation in terms of ordinary non-relativistic 
fluid mechanics [|], ||, |[ ||, 0, ||. The fundamental idea is that sound 
waves propagating in a flowing fluid share many of the formal properties 
of massless scalar fields propagating in a general-relativistic curved space- 
time. Indeed, the propagation of acoustic disturbances in a flowing fluid 
is described by a spacetime metric with Lorentzian signature, the "acous- 
tic metric", which is built up algebraically out of the density, velocity, and 
local speed of sound of the fluid. When the flow is such that there is a sur- 
face where the normal component of the fluid velocity equals the speed of 
sound, the acoustic metric possesses the properties that characterize a black 
hole spacetime in general relativity, and such a surface is therefore called 
"acoustic horizon" . 

As emphasized in [0, 0, acoustic black holes share all the kinematic as- 
pects of relativistic black holes, but do not share in the dynamic aspects. In 
particular, acoustic black holes exhibit Hawking radiation from the acoustic 
horizon, giving rise to a quasi-thermal bath of phonons with temperature pro- 
portional to the "surface gravity" (related to the physical acceleration of the 
fluid as it crosses the acoustic horizon), but they exhibit no simple analog of 
the Bekenstein-Hawking entropy (since that is a dynamical effect intimately 
related to the existence of the Einstein equations in general relativity). 

One of the reasons why acoustic black holes are so popular is that it seems 
that the prospects for experimentally building an acoustic horizon are much 
better than for a general relativistic event horizon. An early estimate can be 
found in [Q] , and related comments are to be found in |^ . Additionally, an 
impressive body of work is due to Volovik and collaborators, who have exten- 
sively studied the experimental prospects for building such a system using 
superfluids such as ^He and ^He g 0, |TT], 0, 0, [TB|, |TB|, [T^, |T§. These 



particular implementations of acoustic geometry make extensive use of the 
two-fluid model of superfluidity, whereas in this paper we will be focussing 
on a conceptually simpler one-fluid model; accordingly, some important tech- 
nical details will differ. For yet another physical implementation of acoustic 
geometries, Garay et al. have investigated the technical requirements for 



implementing an acoustic horizon in Bose-Einstein condensates |19|, and 



some of the perils and pitfalls accompanying acoustic black holes have been 
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discussed in Jacobson's mini-survey pC 



Another attractive feature of acoustic black holes is that they seem to be 
generic, and that they illustrate an important aspect of Lorentz invariance. 
For instance, it is now known, due to the work of Nielsen and collabora- 



tors ||21|, ^ , that in renormalizable non-Lorentz- invariant quantum field 
theories, Lorentz invariance is often an infrared fixed point of the renormaliza- 
tion group equations. Thus, Lorentz invariance can emerge as a symmetry in 
the low-energy limit even if the underlying physics is not explicitly Lorentz in- 
variant. Similarly, in acoustic black holes the underlying physics is explicitly 
classical and Newtonian, but the physics of sound propagation nevertheless 
exhibits a low-frequency approximate Lorentzian symmetry 0. 

In this note we wish to point out a potential difficulty and an opportunity 

— we shall demonstrate that there is a regularity issue that becomes serious 
at the acoustic horizon. Either the Hawking temperature is formally infinite 
(which is the generic situation), or there must be a very precise relationship 
between an external body force that must be applied to the fiuid as it crosses 
the acoustic horizon and the extrinsic geometry of the latter. If this condi- 
tion is not satisfied the "surface gravity" formally diverges, as well as the 
corresponding Hawking temperature. Similarly, the acceleration and density 
gradient of the fiuid at the horizon are formally infinite. For a specified ex- 
ternal force, such divergences are generic, in the sense that they are present 
for almost all fiows, except — in some cases — for a set of measure zero that 
satisfy very special boundary conditions. However, in the case of a constant 
force (including zero force), which is perhaps the most interesting one from 
the point of view of laboratory simulations, no boundary conditions exist 
that correspond to an everywhere regular fiow. 

On the one hand this result suggests that detecting the acoustic Hawking 
effect should be very easy; on the other hand it implies that the naive analysis 
(which demands that both the vorticity and the viscosity be zero) should in 
some way be modified near the acoustic horizon, at least when the external 
forces are such that a formal divergence will certainly occur. For instance, 
adding finite viscosity to the fiuid equations is sufficient in order to regulate 
the surface gravity and Hawking temperature for any choice of external force 

— though finite they can remain large, and can be much larger than naively 
expected. 
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2 Basic equations and assumptions 

The acoustic model of Lorentzian geometry arises from the description of 
the deceptively simple phenomenon of the propagation of sound waves in 
a flowing fluid. Let us therefore recall the fundamental equations of fluid 
dynamics, i.e., the equation of continuity 

^ + f-{pv) = 0, (2.1) 



dt 



and the Euler equation 
where 



pd=f, (2.2) 



dv -> 

a=-^ + {v-V)v (2.3) 



is the fluid acceleration, and / stands for the force density — the sum of all 
forces acting on the fluid per unit volume. We shall assume that the external 
forces present are all gradient-derived (possibly time-dependent) body forces, 
which for simplicity we lump together in a generic term — p V$. In addition 
to the external forces, / contains a contribution from the pressure of the fluid 
and, possibly, a term coming from viscosity. Thus, equation (|2.2| ) takes the 
Navier-Stokes form 




-Vp - p V$ + /viscous, (2.4) 



where 



/viscous = V^^v + (C + ^ ^) V(V • v) (2.5) 

represents the force due to viscous processes, the coefficients rj and ( giving 
the dynamic and bulk viscosity, respectively |^ p5[ . 

In deriving the acoustic geometry, one usually makes a number of tech- 
nical assumptions. 

• The first assumption is that the fiuid has a barotropic equation of state, 
that is, the density p is a function only of the pressure p, so 



p = p{p). 



(2.6) 
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This guarantees that ( |2.1| ) and ( |2.4| ) are a closed set of equations. We 
shall consequentely define the speed of sound as 

^ = ^. (2.7) 

dp 

• The second assumption is that we have a vorticity-free fiow, i.e., that 
V X v =0. This condition is generally fulfilled by the superfiuid 
components of physical superfiuids. 

• A third assumption, often made in the existing literature on acoustic 
geometries, is a viscosity-free fiow. Although this is quite a realistic 
condition for superfiuids we shall see that the presence or absence of 
viscosity can mark a sharp difference in the behaviour of the phonon 
radiation from acoustic horizons. 

These assumptions are sufficient conditions under which an acoustic met- 
ric can be written. However, since the following analysis is independent 
of the introduction of the acoustic geometry (although motivated by it, of 
course), we shall try to be as general as possible, making use of them only 
progressively, as they are needed in order to have an analytically tractable 
system. 



3 Regularity conditions at ergo-surfaces 

Let us start by establishing a useful mathematical identity. If we write v = 
vn , where n is a unit vector and f > 0, then 

dv 

V -v = — + v K, (3.1) 

an 

where d/dn = n • V and K = V ■ n . If the Frobenius condition is satisfied,[| 
then there exist surfaces orthogonal to the fiuid fiow. In this situation, K 
admits a geometrical interpretation as the trace of the extrinsic curvature of 

^ The Frobenius condition is v ■ W x v = 0, or equivalently n ■ \7 x n ~ 0. This 
is sometimes phrased as the statement that the flow has zero "hehcity". The Frobenius 
condition is satisfied whenever there exist a pair of scalar potentials such that v = a V/3, 
in which case the velocity field is orthogonal to the surfaces of constant /3. In view of this 
fact the velocity field is said to be a "surface-orthogonal vector field" . 
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these surfaces. It must be noted that, ahhough zero vorticity is a sufficient 
condition for this to happen, it is not necessary. 

We now focus our attention on the component of the fluid acceleration 
along the flow, an = a ■ n . This can be obtained straightforwardly by pro- 
jecting the Navier-Stokes equation (|2.4| ) along n : 



pan 



2 dp _ d$ - 

P -i- n • /viscous; 

an an 



(3.2) 



where we have used the barotropic condition. 
Next, we rewrite the continuity equation as 



dp dp 



' dv 
dn 



+ vK 



0, 



(3.3) 



noticing that, by the definition (|2.3| ) of a 



where the identity (|3.1| ) has been used. We can express dv/dn in terms of a„ 

(3.4) 



dv 



dv 



at dn 



Thus, equation ( |H73| ) can be rewritten as 



2 dp 2rr , 9v dp 
-V ^-P-K + p--v-. 



p an = —v — pv'K + P — — V — . (3.5) 

Equations ( p. 21) and (|375| ) can be solved for both a„ and dp/dn, obtaining: 



^2 _ y2 



V c K n - J, 

\ dn 



P 



\ 2 ( ^'^ V dp 

VISCOUS I ~l~ C I -;r— — - 



dp 

dn 



^-v'K-^-n.f 

U ±\ lb Jy 

dn p 



\ dv dp 



dt 



dt 



(3.6) 



(3.7) 



In general we see that there is risk of a divergence in the acceleration and 
the density gradient as w — »■ c, which indicates that the ergo-surface^ (the 

^ In general relativity the v —> c surface would be called an "ergosphere" , however 
proving that this surface generically has the topology of a sphere is a result special to 
general relativity which depends critically on the imposition of the Einstein equations. 
In the present fluid dynamics context there is no particular reason to believe that the 
w — > c surface would generically have the topology of a sphere and we prefer the more 
non-committal term "ergo-surface" . 
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boundaries of ergo-regions) must be treated with some dehcacy. The fact that 
gradients diverge in this hmit is the key observation of this paper; we shall 
demonstrate that this has numerous repercussions throughout the physics of 
acoustic black holes. 

Since f ^ = at the ergo-surface it is evident that the acceleration and 
the density gradient both diverge, unless the condition 

d$ 2 1 ^ r dv c dp . . 

c K n ■ /viscous + TTT ^ = (3-8) 

an p ot p ot 

is satisfied on the ergo-surface. Equation ( |3.8| ) is therefore a relationship that 
must be satisfied in order to have a physically acceptable model. Of course, 
it is only a necessary condition, because a„ and dp/dn may diverge at the 
ergo-surface even when ( p.8|) is fulfilled, if the quantities in square brackets 
in the right hand sides of ( |3.6| ) and ( ^.7| ) tend to zero slower than (? — as 
one approaches the ergo-surface. 

For a stationary, non-viscous flow, (|3.8|) reduces to 

^ = cH<, (3.9) 

an 

where again (i<l>/dn, c, and K are evaluated at a generic point on the ergo- 
surface. Thus, in this case it seems that a special fine-tuning of the external 
forces is needed in order to keep the acceleration and density gradient finite at 
the ergo-surface. If the condition (|3.9| ) is not fulfilled but still /viscous = 0, the 
flow cannot be stationary. Near the ergo-surface, an instability will make the 



time derivatives in (|3.8|) different from zero, so that they could compensate 



the mismatch between the two sides of ( |3.9| ). More realistically, we shall see 
later that for a given potential, either no horizon forms, or the flow tries to 
assume a configuration in which (p.9|) is automatically satisfied. 



4 Regularity conditions at horizons 

If we now look at the "surface gravity" of an acoustic black hole it is most 
convenient to first restrict attention to a stationary flow. Defining a notion 
of surface gravity for non- stationary flows is easier in fluid mechanics than 
in general relativity, but is still sufficiently messy to encourage us to make 
this simplifying assumption 0. For additional technical simplicity we shall 
further assume that at the acoustic horizon (the boundary of the trapped 
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region) the fluid flow is normal to the horizon. Under these circumstances 
the technical distinction between an ergo-surface and an acoustic horizon 
vanishes and we can simply define an acoustic horizon by the condition v = 
c. (In complete generality you would have to define an [apparent] acoustic 
horizon as a surface for which the inward normal component of the fluid 
velocity is everywhere equal to the speed of sound 0; this adds extra layers 
of technical complication to the discussion which in the present context we 
have not found to be useful.) Then it can be shown that the surface gravity^ 
qh has two terms P, 0, one coming from acceleration of the fluid, the other 
coming from variations in the local speed of sound. More precisely, qh is 
given by the value attained by the quantity 



'=2 



1 d(c2 



dn 



1 dc^ 

2 dn 



(4.1) 



at the acoustic horizon. (And note that g is defined throughout all space.) 
Now we write, using equations (|373|) and ( p.4|) . 



dc^ 
dn 



dc^ dp 
dp dn 

d^p 
~diP 

d^p 



p K + 



dp2 



Idv" 






V dn^ 


1 V dt 




an 




1 dp 


y2 


dt 


V dt 



(4.2) 



so we find, using (|3.6|): 



^9 p d^p\ / 9 dv V dp 




(4.3) 



at 



Under the present assumptions, time derivatives vanish and v'^ 
the horizon, so it is now evident that the surface gravity (as well as the 
acceleration and the density gradient) diverges unless the condition 



d$ 

dn 



H 



1 

PH 



CrKh n H ■ (/viscous )/f — 



(4.4) 



Hereafter, we label all quantities evaluated at the horizon with the index H. 
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is satisfied. For a non- viscous flow ( |4.4| ) again reduces to ( p.9| ), and the same 
considerations made about the acceleration and density gradient apply. 

Now all this discussion is predicated on the fact that acoustic horizons 
actually form, and would be useless in the case that some obstruction could be 
proven to prevent the fluid from reaching the speed of sound. In order to deal 
with this possibility we shall now check that at least in some speciflc examples 
it is possible to form acoustic horizons under the current hypotheses. For 
analyzing these speciflc cases it is useful to flrst consider generic stationary, 
spherically symmetric flow. 



5 Spherically symmetric stationary flow 

For simplicity, we now deal with the case of a spherically symmetric sta- 
tionary flow in d space dimensions. Spherical symmetry guarantees that the 
fluid flowlines are always perpendicular to the acoustic horizon, and so we 
can ignore the subtleties attendant on the distinction between horizons and 
ergospheres . Additionally, for the time being we shall assume the absence 
of viscosity, /viscous = 0. 

For a spherically symmetric steady inflow, n is minus the radial unit 
vector. Then dn /dn = 0; also 

dn dr' ^ ^ 

and 

K=-^. (5.2) 

r 

Prom equation ( p.3[ ) it follows that a has only the radial component, which 
coincides with — a„ and is 

, c\d - l)/r - d$/dr 

a= -V . (5.3) 

— 

This result could also be obtained directly, without the general treatment of 
section ^. For a steady flow the continuity equation implies 

p V r'^~^ = J = const. (5.4) 



Taking the logarithmic derivative of the above equation one gets 
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On the other hand the Euler equation (|2.4|) takes in this case the form 



p V 



dv 
dr 



dp 
dr 



dr ' 



(5.6) 



where we have used the barotropic condition. Equations ( [5.5| ) and ( p.6| ) can 
be combined to give the useful result 



dv 
dr 



\d-l) 



1 dv 
V dr 



d* 
dr ' 



(5.7) 



which allows one to easily compute the acceleration a = vdv/dr oi the fluid 
for this specific case, recovering equation ( ^.3|) , and to obtain a differential 
equation for the velocity profile v{r): 



dv 
dr 



c2(d- l)/r- (d<l>/dr) 



(5.8) 



When it comes to calculating g, the same analysis as previously developed 
now yields 



g 



2 dp2 / dr 



2 dp2 



r 



(5.9) 



So the acceleration at the acoustic horizon, whose location rn is the solution 
of the equation v{rH)^ = c^rn)'^, formally goes to infinity unless the external 
body force satisfies the condition 



d$ 

dr 



H 



i d -I] 



0. 



(5.10) 



Any further analysis requires one to integrate the differential equation 
(|5.8|). However, this can be done only assigning an equation of state p = p{p), 
and integrating simultaneously equation ( ^.5|) in order to get the dependence 
of c on r. We consider such a specific model in the next section. 



6 Constant speed of sound 

In order to get further insight, let us consider the simple case of a fluid with 
a constant speed of sound. 



Unexpectedly large surface gravities for acoustic horizons? 



12 



It is easy to see that in this case, the condition ( |5.10|) is also sufficient in order 
to keep the physical quantities finite on the horizon. Consider equation ( |5.8|) 
and apply the Bernoulli-de L 'Hospital rule in order to evaluate {dv/dr)^. 
One gets 





+ 



H 



SO {dv/dr)H has a finite value. As a corollary of 
horizon one must have 




d2$ 



dr^ 



< 



(6.2) 

we see that at the 
(6.3) 



H 



H 



SO, in particular, no potential with a non-negative second derivative can lead 
to a horizon on which dv/dr is finite. 

With the assumption ( |6n|) , the differential equation ( ^l8|) for the velocity 
profile can be easily integrated. Its general solution is 



cMn 




2 I 2 

V +Vn 



ro 



(rf-1) In - +$(r) -$(ro), (6.4) 



where ro is arbitrary and vq is the speed of the fluid at ro.0 In order to study 
the general properties of v{r), it is convenient to rewrite equation ( |6.4| ) in 
the form 



where W~'^ is the inverse of the Lambert function |^, deflned as W~^{x) = 
X e^. Given vq and Vq, equation ( |6.5|) implies that the solution v{r) has two 
branches — a subsonic and a supersonic one. This follows immediately from 
the trivial fact that, since W'^^—Vq/c"^) is negative, also W~^{—v^/c^) is 
negative; then, from the plot in flgure ^ we see that there are two possible 
values for v^, one smaller and the other greater than c^. 

We end this section with some remarks that are crucial for a correct 
interpretation of the regularity condition ( ^.10|) . On rewriting ( |6.4| ) or ( |6.5|) 
as 

F{r,v;ro,Vo) = 0, (6.6) 

^Equation ( |6.4| ) simply expresses Bernoulli's theorem. Indeed, it can be written in the 
form + $(r) + h{v, r) = const, where h — J dp/p can be found from ( |5.5| ). 
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Figure 1: The inverse of the Lambert function: W ^(x) = xe^. 

we can represent the location rn of the horizon, for a given potential $ and 
given boundary data (rg, Vq), as the solution of the equation 

F{rH,c]ro,vo) = 0. (6.7) 

On the other hand, differentiating ( |6.6|) and comparing with ( |5.8|) we can 
rewrite the regularity condition ( [5.10D as 

dF 

— {rH,c;ro,Vo) = 0. (6.8) 

It is clear that, if we impose the boundary data (ro,fo), then ( |6.8| ) expresses 
a fine-tuning condition on $ in order to have dv/dr finite at the horizon. 
However, we can reverse the argument and consider the more realistic case in 
which one looks for a physically acceptable flow compatible which an assigned 
$, without trying to force the boundary condition v{rQ) = vq on the velocity 
profile. In this case, equations (|6.71 ) and ( |6l8|) , when solved simultaneously, 
give the location of the horizon, rn, and the value vq of the fluid speed at tq. 
Thus, requiring regularity of the flow for a given potential amounts to solving 
an eigenvalue problem, while if one insists on assigning a boundary condition 
for the speed, a careful fine tuning of $ is needed in order to avoid infinite 
gradients. We stress, however, that although from a strictly mathematical 
point of view both types of problems can be considered, it is the first one 
that is relevant in practice. 
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7 Examples 

We now consider some specific choices, both of $(r) and of v{r), in order to 
illustrate the general situation. 



7.1 Constant body force 

Let us begin with a constant body force, with the linear potential 

$(r) = K r, (7.1) 
where k is a constant. Equation ( |6.4|) becomes, in this case. 



V 

7t. 



c In I — I — f +v^ 



-c' {d - 1) In (-\ + K{r - ro), (7.2) 
Vro/ 



and equation (|6.5|) can be rewritten completely in terms of inverse Lambert 
functions: 



K,r 



2{d~l) 



2{d~l) 



C^{d-1) 



(7.3) 



Following the discussion at the end of section H, we can regard ( |5.10| ) as 
the equation for the locations of th where dv/dr is finite. We have, in this 
case, 

c\d-l) 



(7.4) 



so, excluding the uninteresting possibility th = for d = 1, we see immedi- 
ately that there can be no regular flow with an acoustic horizon when k, < 0. 
For K > 0, one can see that there are no values of vo for which ( [7.4|) is sat- 
isfied. Indeed, setting v = c and r = th = c^{d — 1)/k in ( |7.3| ) gives the 
following equation for Vq: 



1 



-2{d-l) 



(7.5) 



Since, for x < 0, it is > ^(x) > — 1/e (see figure 0), the right hand side 
of equation (|7.5|) turns out to be smaller than —1/e, while the left hand side 
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is greater than — 1/e. Therefore, satisfying equation ( [7.5| ) is impossible, i.e., 
there are no real values of vq that satisfy it, and no regular flow exists in 
which f = c at some point. 

These conclusions are in agreement with equation ( |6.3|) , which implies 
that dv/dr cannot be finite at rfj, because d^$/dr^ = in this case. Thus, 
either v{r) ^ c for all values of r, or dv/dr diverges at the horizon. It is 
not difficult to see that the second possibility is the correct one, because a 
horizon always forms in this type of flow. To this end, let us set f = c in 
(|7.3| ), and look for a solution th (that we do not require to be necessarily 
equal to the one following from the regularity condition, equation ( |7.4| ) — in 
fact, we already know that this would be impossible). We get 



The last factor on the right hand side of this equation is always positive and 
smaller than one, therefore 



W-H-^^] (7.7) 



C2(d- 1 

has the same sign of, and smaller absolute value than, 

(-^) ■ <^-«) 

For K < 0, ( [7.7|) is positive, so also (|7.8| ) is positive and corresponds to a 
positive value rn < tq. For k > 0, ( [TTI ) is negative, so ( |7.8| ) gives two 
positive solutions for th, one smaller and the other greater than tq. In both 
cases, the horizon forms. 

We now illustrate these features by showing some plots of the solution of 
equation (|7.2|) with arbitrarily chosen boundary conditions. Without loss of 
generality we can rescale the unit of distance to set 



K={ 0, (7.9) 




Let us treat these three cases separately. 
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7.1.1 K>0 

For K > figure § clearly shows that there is no obstruction to reaching the 
acoustic horizon. In addition, if we keep the distance scale fixed and instead 
vary k we find the curves of figure ^. 




Figure 2: Plot of the solutions of equation (|6.4|) for several values of d and 
K > 0. We have first fixed k = +c^/ro, and then set c = 1, tq = 1, and 
^0 = 1/2. 

The four things to emphasize here are that: 

1. Velocities equal to the speed of sound are indeed attained; 

2. The gradient dv/dr is indeed infinite at the acoustic horizon; 

3. These particular solutions break down at the acoustic horizon and can- 
not be extended beyond it; 

4. The particular solutions we have obtained all exhibit a double-valued 
behaviour, there is a branch with subsonic flow that speeds up and 
reaches f = c at the acoustic horizon; and there is a second supersonic 
branch, defined on the same spatial region, that slows down and reaches 
f = c at the acoustic horizon. Mathematically, this happens because of 
the double- valuedness of the Lambert function of negative argument, 
as already noted in section ||. 
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Figure 3: Plot of the solutions of equation ( |6.4|) for k = 2, 3, 4 with d = 3, 
c = 1, ro = 1, and vq = 1/2. 

7.1.2 K = (no body force) 

If there is no external body force, then = 1 is uninteresting (the velocity 
is constant). If we now look aX d = 2 and higher then equation (|5.3| ) again 
easily gives us the acceleration of the fluid 



dv n c^id — l)/r 

a = v — = -v^^^ 7-10 

dr — f ^ 



so 

* = _„ ifMzilZr. (7.11) 

dr c — V 

Explicit integration leads us to the solution 

s 1 / 2 2 \ 



which is equivalent to equation (|7.3| ) in the limit k = 0. This can be easily 
plotted for different values of the dimension d as shown in figure ^. 

7.1.3 K<0 

For K < the solutions are plotted in figure ^. 

Finally it is interesting to compare the behaviour of the solutions for the 
different signs of the body force as shown in figure M. 
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Figure 4: Plot of the solution of equation (|6.4|) for k = with the same 
initial values, (c = 1, tq = 1, and Vq = 1/2.) Note the triviality of the d = 1 
solution, which exhibits two branches, with subsonic and supersonic speeds 
respectively. 




Figure 5: Plot of the solution of equation ( |6.4| ) for d = 1,2, 3, and n = 
— c^/ro, where we again set c = 1, ro = 1, and vq = 1/2. 
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In all three cases (k > 0, k = 0, k < 0) we see that the acoustic horizon 
does in fact form as predicted, and that the surface gravity and acceleration 
are indeed infinite at the acoustic horizon. Naturally, this should be viewed as 
evidence that some of the technical assumptions usually made are no longer 
valid as the horizon is approached. In particular in the next section (section 
^ we shall discuss the role of viscosity as a regulator for keeping the surface 
gravity finite. 




0.4 ox 0.8 1 1.2 1.4 IX 1.8 2 2.2 2.4 



Figure 6: Plot of the solutions of equation (|6.4| ) for k 
with c = 1, ro = 1, and Vq = 1/2. 



±1,0 and d = 2, 



7.2 Schwarzschild geometry 

So far, the discussion in this paper has concerned the attainabihty of acoustic 
horizons in general, without focusing on any particular acoustic geometry. A 
more specific, and rather attractive possibility is to attempt to build a flow 
with an acoustic metric that is as close as possible to one of the standard 
black hole metrics of general relativity. Remarkably, this can be done (up to 
a conformal factor) for the Schwarzschild geometry. To be more specific: for 
a fluid with constant speed of sound, one can find a stationary, spherically 
symmetric flow in three spatial dimensions, whose acoustic metric is confor- 
mal to the Painleve-GuUstrand form of the Schwarzschild geometry ^ . This 
possibility has stimulated considerable work concerning the physical realiza- 
tion of an experimental setup that could actually produce such a flow (or, 
more precisely, a two-dimensional version of it |16[)- These particular fluid 
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configurations exhibit a different type of fine-tuning problem tlian tlie one 
we discussed previously. In order to reproduce tfie Painleve-Gullstrand line 



element, tlie speed of tlie fluid must have the profile v = y2M/r, with M 
a positive constant. Then, tq and Vq must satisfy the relation tqVq = 2M, 
and equation ( |6.4| ) allows us to find the external potential needed in order to 
sustain such a flow in d space dimensions: 

«f (r) = «l>(ro) + - + c'(d-l)\n (^) - ^. (7.13) 

Therefore, the potential must be carefully chosen, which will not be easy to do 
in a laboratory. If one does manage to construct such a potential $(r) it will 
automatically fulfill the fine-tuning condition (|3.9| ) at the acoustic horizon. 



rn = 2M/c^. This is only to be expected, because df/dr = —^M/2r^ blows 
up only as r ^ 0. Also, since we know the surface gravity of a Schwarzschild 
black hole is finite, any fluid flow that reproduces the Schwarzschild geometry 
must by definition satisfy the fine tuning condition for a finite surface gravity. 
Looking at the issue from the point of view discussed at the end of section 



one expects that, given the potential ( |7.13| ), the value Vq = J2M/ro 



is the solution of equations (|6.7| ) and ( |6.8|) , while f(r) = y2M/r is the 
corresponding eigenfunction that is selected by the requirement of having a 
regular flow. This is indeed the case: Equation (|6.7|) now gives = 2M/ 
which, substituted into (|6.8|), leads to the following equation for Vq: 



2 \roV^ 




It is trivial to check that vq = J2M/ro is, in fact, a solution of ( [7. 14|) . 



Considering the same potential (|7.13|) , but values of Vq different from 



2M/ro, corresponds to flows either with no horizon, or in which {dv/dr)^ 
diverges. This is evident in figure 0, which confirms the "eigenvalue char- 
acter" of the problem of finding a regular flow. Notice that there are two 
solutions that are regular at the horizon, with opposite values of {dv/dr)^, 
in full agreement with the fact that equation ( |6.2| ) only determines the square 
of {dv/dr)H- 

Additionally, note that what we have done above has been to ask how to 
mimic a slice of the (3 + l)-dimensional Schwarzschild geometry with a (rf+l)- 
dimensional fluid flow. We could ask what happens in different spacetime 
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Figure 7: Plot of the solutions of equation (|6.4| ) for the potential ( |7.13|) , 
with M = 1/2, d = 3, a = 1, and ro = 4. The corresponding boundary 
conditions are Vq = 0.4, vq = 0.5, and Vq = 0.6. 

dimensions: for the {D + 1) -dimensional generalization of the Schwarzschild 
geometry the fluid flow generalizes to v = ^j2M/r^~'^, and the potential 
gradient required to produce this flow is 

$(r) = $(ro) + ^2+^\d--)\n (-] - (7.15) 

Again a very specific external body force is needed to set up the very specific 
fluid flow corresponding to a higher- dimensional Schwarzschild geometry. 

7.3 Reissner— Nordstrom geometry 

We mention in passing that generalizing this discussion to the (3 + 1)- 
dimensional Reissner-Nordstrom geometry is straightforward. This geom- 
etry is described in Painleve-GuUstrand form by the fluid flow 




(7.16) 



The external potential required to set up this fluid flow is then 

$(r) = (d - — + + - + const. (7.17) 

^ ' V 2/ Vro/ r 2r2 2 V 2Mr ) ^ ' 
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7.4 The canonical acoustic black hole 

To wrap up our section on specific examples, we add a few words about the 
"canonical" acoustic black hole discussed in ^. In that model the fluid 
is assumed to have a constant density throughout space, and the continuity 
equation is then used to deduce the velocity proflle 

v{r) oc l/r^. (7.18) 

Note that "constant density" is actually a much weaker statement than in- 
compressibility, and the word incompressible should be excised from all of 
section 8 of and replaced by this phrase. Now in 0, || the velocity proflle 
was determined purely on these kinematic grounds, and no attempt was made 
to put this background fluid flow back into the Euler equations to determine 
the external body force required to set up the flow. (In that paper, almost all 
the attention was focussed on the fluctuations rather than the background 
flow.) 

Determining the potential is an easy application of the general analysis 
of this article [see equation (|5.71 )]. We calculate 

$(r) oc + const. (7.19) 

With hindsight this can be seen to be nothing more than a special case of 
Bernoulli's theorem for a constant-density flow 

$(r) = -1 + const. (7.20) 

The single over-riding message coming from all these speciflc examples is 
the generic dichotomy between a formally inflnite surface gravity and need- 
ing a highly speciflc boundary condition to be satisfled. In the following 
section we shall regulate the generically inflnite surface gravity by using a 
less idealized model for the fluid. 

8 Viscosity 

A viscous flow is governed by the Navier-Stokes equation ( |2.4| ). In general, 
there will be two contributions to viscosity, associated with the coefficients 
rj and (. Since our treatment in the present section does not pretend to be 
realistic, but we simply wish to point out how viscosity acts as a regulator 
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for tlie surface gravity, we shall set the bulk viscosity coefficient C to zero, 
in order to have a model with as few free parameters as possible. (This is 
sometimes called the "Stokes assumption" ||2^.) For a spherically symmetric 
inflow one has 



n 



V'v +-W -v] = - {— + 



4 
3 



dr 



.11 



so the Navier-Stokes equation 
dv 



becomes, in the stationary case. 



dr 



^2 _ y2 



c\d-r 



d<^ 

dr 



+ 



Au (d^v {d-l)dv {d-l)v" 



dr^ 



+ 



dr 



(8.2) 



where we have introduced the coefficient of kinematic viscosity v = rj/p. 
Hereafter we assume, for simplicity, that is a constant. (However, any hy- 
pothesis about u must ultimately by justified by a kinetic model for the fluid, 
and it is worth noticing that there are plausible distribution functions that 
lead to a velocity-dependent u; see, e.g., [^. Obviously, such a dependence 
can have important repercussions on the conclusions of the present section.) 
The acceleration is then infinite at the horizon unless 



Au \ c^{d - 1) 



3cr 



H 



d$ 
dr 




+ 



{d 



H 



dv 
dr 



0. (8.3) 



Since this involves higher-order derivatives at the horizon, it can no longer be 
regarded as a fine-tuning constraint, or as an equation for r/^, but merely as 
a statement about the shape of the velocity profile near the horizon. Indeed 
the general solution to this equation when c is constant is 



v{r) 



c + 




r - rn 



H 




-aH\[r - rn 



+ 



with an arbitrary, and finite, at the horizon. We can rearrange 
a differential equation for v{r): 



(8.4) 
to get 



fv 
dr'' 



+ 



■(d-i: 



{d-l)v 



4z/ 




c\d-iy 



15) 
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Unfortunately, integrating this equation is completely impractical in general 
and we must resort to the analysis of special cases. 

8.1 d = 1, constant body force 

Even the case of constant body force is intractable unless d = 1, in which 
case we get (following the steps above) 



df 
dr 



-K + 



4z/ d'^v' 
3 dr^ 



(8.6) 



This single second-order differential equation can be turned into an au- 
tonomous system of first-order equations 

df 
dr 



^ = n, 



dn 

dr 



3k 3v 



17) 



1 - 



n. 



We can plot the fiow of this autonomous system in the usual way and it 
clearly shows that it is possible to cross the acoustic horizon f = c at arbitrary 
accelerations an and arbitrary surface gravity qh (see figure 



8.2 d = 1, zero body force 

Integrating equation (EM) once (this is easy provided c is a constant), we get 



df 
dr 



df 
dr 



2 I 2 

f +fo 



+ -(<l>(r)-$(ro)), (8.8) 



where (ro,fo) again denotes an arbitrary pair of initial values. If d = 1 and 
K = this equation for df (r)/dr reduces to 



df 
dr 



df 
dr 



3 

8^ 



2 I 2 

f +fo 



(8.9) 



In this particular case the analysis is sufficiently simple that we can say 
something about the acceleration at the horizon, namely 



df 
dr 



H 



df 
dr 



^.10) 
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Figure 8: Plot of the phase space for d = 1, k = constant > 0, and nonzero 
viscosity. The transverse hne identifies a separatrix in the integral curves. 
Note that the integral curves can intersect the acoustic horizon at arbitrary 
values of the surface gravity. 



That is 



dv 

Qh = an = c — 
dr 



'■o 



3c 



2 I 2 
C +Vq 



(8.11) 



which is an explicit analytic verification that viscosity regularizes the surface 
gravity of the acoustic horizon. 

We can plot the flow in the usual way and it again clearly shows that it 
is possible to cross the acoustic horizon = c at arbitrary accelerations an 
(see figure ^. 



8.3 (i > 1, constant body force 

The relevant equation is 



dr^ 



'(d- 1) ^ 3 -v'^\ dv {d-l)v 



c2(d- 1) 



4z/ 



dr 



4z/ 



(8.12) 
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Figure 9: Plot of the phase space for d = 1, k = 0, and nonzero viscosity. 
Again note that the integral curves can intersect the acoustic horizon at 
arbitrary values of the surface gravity. 

which can be recast as 

f ? = n 

dr 

dn ^ 3 / c\d-l) \ ^ {d-l)v ^ f?>vf 

dr Av \ r J yAu \ v"^ J r J 

(8.13) 

This is no longer an autonomous system of differential equations, (there is 
now an explicit r dependence in these equations) so a flow diagram is mean- 
ingless. Nonetheless the system can be treated numerically and curves plot- 
ted as a function of initial conditions. As an example we plot some curves 
in the phase space ioi d = 2, and verify that at least some of these curves 
imply formation of an acoustic horizon. 

As a final remark we think it is useful to briefly discuss the effect of 
viscosity with regard to Hawking radiation. It has been shown that the 
addition of viscosity to the fluid dynamical equations is equivalent to the 
introduction of an explicit violation of "acoustic Lorentz invariance" at short 
scales P, 0. Thus one may wonder if such an explicit breakdown would 
not lead to a suppression of the Hawking flux as well. Indeed the violation 
of Lorentz invariance is important for wavelengths of order Aq = z^/c ||^, 
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Figure 10: Plot of some solutions to the non-autonomous system (R7 



for various initial conditions for d = 2, k = constant, and nonzero viscosity. 
Note that at least some of these curves intersect the acoustic horizon, and 
do so at various but finite values of the surface gravity. 



0, introducing in this way a sort of cutoff on short wavelengths which can 
dramatically affect the Hawking flux 0, |^, |^, |2y] . 

Thus there is naively a risk that using viscosity to remove the unphys- 
ical divergences at acoustic horizons would also "kill" the phenomenon one 
is seeking. This problem has been extensively discussed in the literature 
(see e.g. [|^) and it has (quite remarkably) been demonstrated that such 
a violation of Lorentz invariance is not only harmless but even natural and 
useful. In particular, viscosity can be shown to induce P, |^ the same type of 
modifications of the phonon dispersion relation which are actually required 
for circumventing the above cited problem [|], |], ^ So the emergence of 
viscosity appears to be indeed a crucial factor, both for allowing the formation 
of acoustic horizons and, at the same time, for implementing that mechanism 
of "mode regeneration" which permits Hawking radiation in presence of short 
distance cutoff. 
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9 Conclusions: Danger+opportunity. 

Let us summarize the results that we have obtained for a fluid subjected 
to a given external potential: In the viscosity-free, stationary case, we have 
seen that if the flow possesses an acoustic horizon, the gradients of physical 
quantities, as well as the surface gravity and the corresponding Hawking flux, 
generically exhibit formal divergences. There are two ways in which a real 
fluid can circumvent this physically unpalatable result. For a broad class 
of potentials, there is one particular flow which is regular everywhere, even 
at the horizon. In this case, it is obvious that the fluid itself will "choose" 
such a configuration. Mathematically, imposing the regularity condition at 
the horizon amounts to formulating an eigenvalue problem. However, there 
are physically interesting potentials — such as a linear one — for which this 
is impossible. We view this result as both a danger and an opportunity. 
A danger because infinite accelerations are clearly unphysical and indicate 
that the idealization of considering a irrotational barotropic inviscid perfect 
fluid (and this idealization underlies the standard derivations of the notion 
of acoustic metric [|I], ||, ^, |[), is sure to break down in the neighborhood of 
any putative acoustic horizon. Indeed, the divergences will be avoided in real 
life simply because one or more of the simplifying hypothesis become invalid 
as f — i> c. 

On the other hand, this may be viewed as an opportunity: Once we 
regulate the infinite surface gravity, by adding for instance a finite viscosity, 
we find that the surface gravity becomes an extra free parameter, divorced 
from naive estimates based on the geometry of the fluid flow. The common 
naive estimates of the surface gravity take the form ^ 

9H ^ ^, (9.1) 

where i? is a typical length scale associated with the flow (a nozzle radius, 
or the radius of curvature of the horizon). The analysis of this note suggests 
that this estimate may in general be misleading because it does not take into 
account information regarding the dynamics of the flow. Because of this the 
surface gravity could be considerably larger than previously expected. 

There is a potential source of confusion which we should clarify before 
wrapping up: in general relativity the physical acceleration of a stationary 
observer hovering just outside the black hole horizon diverges, but when 
an appropriate red-shift factor is applied and the properly defined surface 
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gravity is calculated that surface gravity proves to be finite. On the contrary, 
in the acoustic black holes it is the Newtonian acceleration of the infalling 
observers that (in the absence of fine-tuning) diverges at the horizon, leading 
to an infinite surface gravity. Why the difference? It is here that the actual 
dynamical equations governing the background geometry come into play. 
The physics that is identical between gravitational black holes and acoustic 
black holes is the kinematical physics of fields propagating in the respective 
Lorentzian spacetime metrics. The physics which is different is that which 
depends on the dynamical equations of motion of the background geometry. 
For gravity, the latter is governed by the Einstein equations while for acoustic 
black holes it is governed by the hydrodynamic equations — these equations 
are sufficiently different that the geometries of the two Lorentzian metrics 
can be quite different, even though qualitative features such as the existence 
of event horizons may be quite similar. 

Our general discussion, plus the specific example utilizing viscosity, makes 
it clear that it is the specific technical restrictions placed on the hydrody- 
namic equations that lead to the formally infinite surface gravity — and so 
one might wonder how much of the current analysis to trust. For exam- 
ple, in real superfluids the existence of roton excitations leads to a break- 
down of irrotational flow before the acoustic horizon is reached . Adding 
vorticity is certainly technically complicated (see for instance the recent 
book by Ostashev p8[), but this may merely be a technical complication. 



not a fundamental barrier to progress. For technical discussions regard- 
ing the possibility (probability) of actually building acoustic black holes 
see i, 0, II, llll 111, |18|, PI, ll. Note that the Garay et 
al. implementation of acoustic black holes |19[ is built on a different physi- 
cal background; they use Bose-Einstein condensate governed by the Gross- 
Pitaevski equation rather than a barotropic fluid governed by the Euler- 
continuity equations. Therefore the perils and opportunities delineated in 
this article do not necessarily apply to their particular situation. A simi- 
lar remark applies to Volovik's implementation based on two-fluid models 
of superfluidity (for example ^He-A), where the horizon is defined using the 
speed of the quasi-particles, rather than by the speed of sound per se. Of 
these two speeds, the former is much smaller than the latter, so the surface 



gravity at such horizons is always finite [18|. In short, while specific physical 
implementations of the acoustic geometry idea all have their characteristic 
peculiarities and potential pitfalls, overall the experimental prospects con- 
tinue to look extremely promising. 
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